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ABSTRACT 

In the case of an initially conical jet, we study the relation between jet collimation by the 
external pressure and large-scale morphology. We first consider the important length-scales 
in the problem, and then carry out axisymmetric hydrodynamic simulations that include, for 
certain parameters, all these length- scales. We find three important scales related to the col- 
limation region: (i) where the sideways ram-pressure equals the external pressure, (ii) where 
the jet density equals the ambient density, and (iii) where the forward ram-pressure falls below 
the ambient pressure. These scales are set by the external Mach-number and opening angle of 
the jet. We demonstrate that the relative magnitudes of these scales determine the collimation, 
Mach-number, density and morphology of the large scale jet. Based on analysis of the shock 
structure, we reproduce successfully the morphology of Fanaroff-Riley (FR) class I and II 
radio sources. Within the framework of the model, an FR I radio source must have a large 
intrinsic opening angle. Entrainment of ambient gas might also be important. We also show 
that all FR I sources with radio lobes or similar features must have had an earlier FR II phase. 
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Extragalactic radio sourc es are traditionally divided into two 
morphological classes (iFanaroff & Riley 1 9741) : those in class I 

' (FR I) are edge- d arkened and those in class II edge-brightened. 

' iFanaroff & RilevI found that the two classes differ in radio power, 
those with 178 MHz luminosities below 2 x 10^^ W Hz"^ sr~^ 
(Hubble's constant = 50 km s~^ Mpc~^) being in class I and those 
above in class IL In subsequent studies, the critical power, Pcrit^ 
was found to be corr elated with the optica l magnitude of the host 
galaxy {P^^ qc lI^,, |Ow en & Laing' "1989"; 'Owen & Whit elll99ll : 
IOwenill993l : [Ledlow & Owen 1996; Ghisellini & Celotti 2Q0i|) It 
is well known that brighter gal axies tend to have bi gger central 
super-massive black holes (e.g. iGiiltekin et al.l l2009b . which are 
responsib le for the jet production in th e first place. Interestingly 
however, IWold, Lacv, & Arrnusl (l2007h do not find a correlation 
of FR class with the mass of the super-massive black hole in the 
host galaxy for a sample of 21 3CRR radio sources near the criti- 
cal power. Thus, in determining the large scale morphology, the jet 
power is the most important influence - at least one other also plays 
a role, but this factor does not seem to be the black hole mass. 

FR I radio sources often show narrow jets emanating from 
the core, which is usually thought to coincide with the AGN, and 
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extending out to a bright flaring point. Beyond this flaring point, 
the flow widens and dims. These findings have been modeled 
as a t ransitio n from a laminar to a tur bulent flow (e.^. i Bicknelll 
11984 Il986al lbl: [K^issarovlll990 dlbllal: IWang et al.ll2009l) . which 
mixes turbulently with the entrained ga s. This basic mod el has 
been confirmed by numerical simulations: iRossi et al.l(l2008h show 
that a jet that is forced t o entrain may make the desired transition. 
IPerucho & Martil (l2007h find the transition in flow character after 
an induced strong expansion phase in the jet. The reason for such 
an expansion phase and why it should occur at low jet power only , 
rem ains however unc lear. iHardcastle. Sakelliou, & Worralll (l2005h 
and lJetha et al.l (l2005h present studies of morphological features in 
three FR I sources, focusing on the density structure in the ambient 
X-ray gas as observed by Chandra and XMM-Newton. They find 
the flaring point is located at a place where the ambient gas temper- 
ature rises steeply from the inner cool core to the hotter intra-cluster 
gas. This is accompanied by a local flattening of the density profile. 
In each case, the location is a few tens of kpc from the core. It is 
clear that observations suggest the FR class of a radio source is 
significantly influenced by both the intrinsic jet properties and the 
properties of the environment. 

The physics of the structure and evolution of FR II sources 
has been studied in some detail. Such sources often show evi- 
dence of a jet at various places from the core out to the edge 
of the source, which might be over a Megaparsec away (e.g. 



2 M. Krause et al. 



iMullin. Hardcastle. & Rilevl2006h . There, they feature one or more 
hotspots. A diffuse, cyHndrical "cocoon" extends from the hot 
spots back towards the radio core. Axis ratios of the se cocoons 
vary between one an d about six (e.g. iKaiser & Alexander 1997i. 
and references therein. iMulHn. Rilev, & Hardcastlel20"08l) . Th is has 
inspired self-similar models of FR II radio source growth (iFalld 
ll99ll : lKaiser & Alexandeilll997h . These models are valid far away 
from the characteristic sca les of the problem. An outer scale L2 
(iKomissarov & Fallell998h . at around which the source comes into 
pressure equilibrium with its surroundings, is given by 
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where Qq is the kinetic power of the source, px is the external 
density, and Cx is the external sound speed. F R II radio sources 
approaching L2 have been studied analyticallv (lAlexandei l2002h 
and n umerically fe.g. lKraus3l2005l : rGaibler, Krause. & Camenzindl 
l2009h . Up to L2, the source is overpressured with respect to its en- 
vironment, and is also expected to have a strong bow shock (weak- 
ening as the source approaches L2). Observati ons of weak bow 
shocks in radio sources ^ 100 kpc in size (e.g. iMcNamara et al.l 
I2OO5I : iNulsen etal.ll2005h may be simply understood in terms of 
them being close to L2. 

Crucial to the existence of the self- similar solution is the self- 
confinement of the jet by the source's own cocoon pressure. This 
happens automatically in numerical models of the propagation of 
light jets (heavy jets show only a rudimentary cocoon). In such sim- 
ulations, the overpressured cocoon drives a re-collimation shock 
into the beam as soon as the latter enters the computational domain, 
rega rdless of whether the beam is initially conical or cylindri cal 
(e.g. IKomissarov & Fallel[T998l : lKrausd[2003l : lHeath et"ai]l2007h . 

An initially overdense, conical beam, may be collimated via 
a re-collimation shock. It occurs when the external pressure px be- 
comes comparable to the sideways ram-pressure, pjv? sin^ 6, where 
Pj is the jet density (oc in a conical jet, where r is the distance 
to the AGN), v, th e cons tant beam velocity, and 6 the half opening 
angle JScheuei 's (1 19741) model B). This position is rel ated to and 
occur s somewhat downstream of the inner scale Li fe.g. lAlexandeij 
l2006h . given by 
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On scales Li <C L <C L2, we expect self- similar behaviour. At Li 
the jet density becomes roughly comparable to the external density, 
and the mass of the swept-up gas becomes comparable to the mass 
that has gone through the jet channel. 

The development of the large-scale morphology of the ra- 
dio source depends upon what happens around Li. If a proper re- 
collimation shock is driven into the beam, an FR II source will de- 
velop. Otherwise, a morphology of the FR I type might result. 

Here, we present jet simulations on a spherical grid starting 
at a fraction of L] out to a fe w hundred times Li. We follow the 
approach of (lAlexandei]|2006h in employing non-relativistic hydro- 
dynamics, addressing the same scale, downstream of Li. This al- 



lows detailed studies of the re-collimation process and its relation to 
the large-scale morphology. We show that a proper re-collimation 
shock only occurs for small opening angles. 

We further discuss the important scales of the problem in Sec- 
tion [3 Numerics and simulation setups are presented in Section [3] 
The simulation results are shown in Section |4] and discussed in 
Section (5] We conclude in Section [6l 



2 RE-COLLIMATION, COCOON FORMATION AND 
TERMINAL SHOCK LIMIT 

The physics of the region downstre am of L] involves understand- 
ing of the three length- scales (see lAlexanded (l2006h for further 
discussion. Table [T] summarises the definitions), associated to re- 
collimation (Lia), cocoon formation (Lib) ^nd the terminal shock 
(Lie). The first, Lia, is where the jet's sideways ram-pressure equals 
the ambient pressure leading to potential re-collimation by the am- 
bient gas: 
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where 7 is the ratio of specific heats in the environment, Q. = 
l7i(l — COS 0), the solid angle of the jet, and the Mach-number 
with respect to the ambient sound speed. If, however, an overpres- 
sured cocoon forms, the jet-collimation can occur on smaller scales. 
Before the jet collimates, the jet density decreases with distance due 
to its fixed opening angle. The jet starts to develop a cocoon at the 
second length-scale Lit, where the density of the jet falls below that 
of its environment: 
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The third. Lie is the length-scale where the ram-pressure of the 
un-collimated conical jet would fall below the pressure in the envi- 
ronment (see below for the case when an overpressured cocoon is 
present): 
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Since the terminal shock is advanced by the jet's forward ram- 
pressure, it will advance only up to Lie, unless the jet is colli- 
mated before reaching this third scale. If an overpressured cocoon 
is present, the critical scale for the jet termination shock will be 
further in. However, the scales for re-collimation and jet termina- 
tion will be affected in the same way (compare eqs (3) and (5), the 
dependence on the environmental pressure is via the sound speed 
which enters the Mach number and for this consideration the en- 
vironment would be the cocoon). The terminal shock is the prime 
site of particle acceleration and therefore synchrotron emission. It 
may therefore be identified with the hotspots of a radio source. If 
it remains at Lie because the jet fails to collimate beforehand, the 
cocoon may still expand and the stationary shock at Lie n^^Y be 
identified with the flaring point of an FR I source. 

It is interesting to look at the relation of these scales to one 
another. The ratio Li^/Lia is given by: 
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If Lla < ^ib then the jet will be collimated by the external pressure 
before the jet density has fallen to that of the surroundings. Once 
collimated the jet density does not decline further along its length. 
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Table 1. Definition of length-scales 



Length-scale fonnula symbol assoc. transition^ 

/ \ 1/2 

Inner (^j L. 

Re-coUimation y^/^Mx sin0Li/(2^2^/2) 

Cocoon formation L\/{2Q}/^) Lib -^ib 

Terminal shock y^l'^M^Li / {IQ}/'^) Lie -^ic 



Outer ( ^ ) ' L2 



Symbol of the associated transition point when the jet is in a cocoon. 



Therefore, for Lia < Lib ^ cocoon does not form. Conversely, for 
Lia > Lib ^ cocoon is formed first and the jet is collimated by the 
pressure in the cocoon. Th is s econd case is stud ied in detail by 
iKomissarov & Falld (Il998h and lAlexandel (l2006h . and the source 
transitions to near self- similar subsequent evolution. The ratio 

gives the ratio of the collimated jet density to the external density. 
For a cocoon to form, T] < 1, and for a substantial cocoon, T] << 1. 
From Very Long Baseline Interfer ometry (VLB I), it is clear that 
the jets are initially relativistic fe.g. iBritzen et al.lEoOSh : on scales 
of ^ 100 kpc, their velocities are still a reasonable fraction of 
the sp eed of light fe.g. lHardcastle et al .111 9991 : iMullin & Hardcastll 
\200% . Hence, the external Mach-number for a jet is typically of 
order several hundred. For very small opening angles of 1 de- 
gree, therefore a cocoon will not form. However, for larger opening 
angles cocoons will form. 

The ratio Lic/Lib is given by: 

W^lb = V^Mx. (8) 

If Lic/Lib >> 1 then a lobed source, i.e. one that has formed a 
cocoon, also has the hotspots near the leading edges. From Equa- 
tion ([H), this is always the case. Hence, every radio source that has 
lobes went through a phase when it had its hotspots near the leading 
edges, i.e. it was once an FR II source. To remain an FR II source, 
the re-collimation shock needs to reach the axis, before th e je t run s 
out of thrust. Calculation of this point is difficult. Alexander (2006) 
argues that re-collimation, in the case when the radio source is col- 
limated by the pressure of a cocoon, should be completed at two to 
three times a length scale Lr. Lr corresponds to Lia in the present 
paper, and the numerical factor is in excellent agreement with the 
results of the simulations we present below. Adopting a value of 
2.5, the critical ratio becomes: 

Lic/2.5Lia = 1/(2.5 sine). (9) 

Thus, jets with a half opening angle below about 24° have enough 
thrust to propel the terminal shock past the re-collimation shock. 
Once this occurs, the terminal shock may continue essentially in- 
definitely, since the beam is collimated, and the ram-pressure stays 
constant (the exception being an inverted pressure profile in the ex- 
ternal medium). Such sources will end up as large-scale FR II radio 
sources. 

Jets with a half opening angle greater than 24° have the re- 
collimation shock so far downstream that the terminal shock is un- 
able to reach it. In this case, the hotspot will stall at Lie, and might 
turn into a flaring point. The exhaust plasma that leaves the beam 



may still form a classical cocoon for some time. The plasma down- 
stream of the flaring point streams trans sonically towards the tip of 
the cocoon. There is no clear distinction between jet and cocoon 
in this region. Ram-pressure balance between this plasma and the 
shocked ambient gas gives the velocity Vcd of the contact surface 
as Vcd = \/W^coc where T]' is the density ratio between the cocoon 
and the shocked ambient gas, and Vcoc is the bulk velocity of the 
cocoon plasma. Therefore, as long as the cocoon plasma is under- 
dense with respect to the shocked ambient gas, the contact surface 
must advance more slowly than the cocoon plasma. The cocoon 
plasma therefore builds up thermal pressure near its leading edge 
and forms a backflow in a part of the cocoon. The situation changes 
at some point because of entrainment of shocked ambient gas into 
the cocoon. When enough dense gas is mixed into the cocoon so 
that - spatially averaged - T]' ^ 1, the contact surface advances at 
the same velocity as the cocoon plasma. The morphology is then 
transformed into a chimney-like outflowing structure. 



3 NUMERICS 

We have carried out 2.5D axisy mmetric hydrodynamics simula- 
tions with the FLASH code (v. 2.5 lFrvxell et al.l2000l : ICalder et all 
l2002h to demonstrate the physics outlined in the previous sec- 
tion. FLASH evolves the equations of compressible hydrodynam- 
ics conservatively, using third order (piecewise parabolic) interpo- 
lations. We have run six simulations of conical jets into a constant 
ambient density, chosing units of Li (length), px (external density), 
and Cx (external sound speed). The ratio of specific heats is 7= 5/3. 
In these units, the external pressure is therefore 3/5. 

We use spherical polar coordinates and the jet is implemented 
as a boundary condition. Other boundaries are reflective. In Sec- 
tion [21 we found that the opening angle and the Mach-number de- 
termine the characteristics of the flow. Therefore, we vary these pa- 
rameters in the simulations. The other parameter that is required to 
define the simulation setup is the pressure ratio between the jet and 
the ambient medium. We set it to unity at the jet injection bound- 
ary, noting that the adiabatic expansion renders the beam pressure 
quickly insignificant. Finally, the computational domain specifies 
the outcome of the simulation. Our aim is to include all the relevant 
scales of Section [2l However, choosing a realistic external Mach 
number of 500, and starting the jet smaller than Li , we were unable 
to run the simulation for long enough for the jet to reach L2, as three 
to four orders of magnitude between the smallest and largest scales 
are required for the three opening angles we chose. Making use of 
the polar coordinates and the adaptive mesh refinement, we reach a 
few hundred times Li . Therefore, we have chosen a second set of 
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Table 2. Simulation parameters 



Run 


















•^Ici/^l 




dr^U 




M05 


500 


5 


0.363 


210 


182 


11 


3.23 


2087 


130 


3953 


0.142 


0.088 


M15 


500 


15 


0.125 


76 


181 


11 


1.08 


697 


41 


3953 


0.049 


0.352 


M30 


500 


30 


0.062 


50 


176 


9 


0.54 


352 


18 


3953 


0.024 


0.352 


m05 


5 


5 


0.262 


22 


1.82 




3.23 


20.87 




3.95 


0.030 


0.044 


ml5 


5 


15 


0.098 


18 


1.81 




1.08 


6.98 




3.95 


0.039 


0.022 


m30 


5 


30 


0.050 


16 


1.76 




0.54 


3.52 




3.95 


0.019 


0.059 



^ External Mach-number 

^ Half opening angle / degree 

^ Inner boundary of the computational domain 

^ Final size 

^ Length-scale where the sideways ram pressure equals the ambient pressure (re-coUimation) 

^ Length-scale where the sideways ram pressure equals the cocoon pressure at the end of the simulations 

^ Length-scale where the density of the uncoUimated jet would equal the ambient density (cocoon formation 

^ Length-scale where the forward ram pressure of the uncoUimated jet would equal the ambient pressure (terminal shock limit) 

^ Length-scale where the forward ram pressure of the uncoUimated jet would equal the cocoon pressure at the end of the simulations 

^ Outer scale 

^ Radial resolution (max.) 

^ Meridional resolution (max.) / degree 



M05, time =2.26 i ^ m05, time =6.46 ^» 




Figure 1. Final logarithmic density (top half) and pressure (bottom half) distributions for all simulations. The time unit is L\ divided by the ambient sound 
speed. The plots use different x and _y scales in order to make the beam features visible. For run m30, the density scale saturates at ten to enhance contrast, 
which affects the beam near the injection region, only. 



simulations with an external Mach-number of 5. For this second set, 
all the relevant scales, including L2, are within the computational 
domain, and we reach four to five times L^. We choose three open- 
ing angles around the critical one defined in Equation (|9j. We label 
the simulations with a capital (small) M for the high and low Mach- 



number cases, respectively, followed by the half opening angle in 
degrees. The simulation parameters are summarised in Tabled 
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4 RESULTS 
4.1 Structure 

The spatial structure of the simulated radio sources is shown in 
Figure [T] From the discussion above, and Table |2] we expect - and 
find - all the high Mach-number simulations to be cocoon domi- 
nated, since L^, is only a few Li, and Lia, the length-scale which 
determines collimation, is much greater. L2 is 20 to 80 times the 
maximum jet size, and therefore all the three jets are overpressured 
with respect to the external medium throughout the simulations. 
The strong bow shock seen in all three cases is evidence for this 
(Figures |2]and[3]). Because of the high cocoon pressure, collimation 
occurs before Li^\ there is however reasonable agreement with the 
length scales calculated for the sideways ram-pressure equalling 
the cocoon pressure. We denote the scales which refer to the co- 
coon pressure at the end of the simulations by xia and xic. They are 
also given in Table|2] Since the cocoon pressure declines with time, 
the re-collimation shock reaches the axis later for later simulation 
times (Figure El X ^ lOOLi at t = 1.35, X ^ 140Li at t = 2.71 
and X ^ I8OL1 ait = 5.42, note the varying scale on the X-axis ). 
The re-collimation shock still reaches the axis comparatively late 
in run M05 - this is due to a pressure gradient in the cocoon to- 
wards the left (Figurejl}. Th ese features have already been noted by 
iKomissarov & Fall^ (Il998h . and we find good agreement with their 
simulations. As predicted by Equation (|9]), the re-collimation shock 
always reaches the axis for the simulated jets with small opening 
angles and never reaches it for the ones with the large opening an- 
gles. 

All the high Mach-number runs form very light jets and co- 
coons with T] ^ 10~^. This is consistent with their Li^Aia ratios. 
At the times shown in Figure[T] the cocoons are overpressured by a 



factor of several tens to hundreds. As the jets evolve towards L2, the 
cocoon pressure drops and the collimated jets become even more 
underdense. 

As predicted in Section|2] for run M30 the terminal shock can- 
not be propagated through the re-collimation shock. The jet does 
not therefore re-collimate properly, and a stable collimated beam 
cannot form. The terminal shock is however not strictly a standing 
shock. Its location is where the forward ram pressure equals the 
cocoon pressure (xic). At the end of the simulation xic = I8L1, as 
calculated from the current cocoon pressure. This is in good agree- 
ment with the measured position. As can be seen from Figure |3] 
the position of the terminal shock advances with time, according 
to the declining cocoon pressure. In our simulation, the terminal 
shock never advances beyond the re-collimation shock, as expected 
from Equation ^ - the ratio between the terminal shock limit and 
the re-collimation scale only depends on the initial opening angle; 
run M30 demonstrates this result for a wide range of cocoon pres- 
sures. 

The low Mach-number runs may be found in the right column 
of Figure [1] In run m05, the re-collimation shock reaches the axis 
at about 3 Li, consistent with its Lia value. This is only slightly 
larger than Lib, ^nd hence essentially no cocoon is formed. We 
expect - and find - the first reflection point of the re-collimation 
shock, where it meets the axis, to have a constant x- value through- 
out the simulation (Figure |4]). The density structure of the beam is 
a hollow spine with a dense sheath. On average, the beam stays 
overdense with respect to its environment, as expected from Equa- 
tion ([7]). A minor filamentary cocoon forms for run ml 5 (Figure (TJ. 
The beam density in this case is almost matched to the environ- 
ment. This is also as expected from the values of Lia and Lit, (Ta- 
ble |2]). For run m30, Li^/Lia = 0.3, and we expect an underdense 
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Figure 4. Logarithmic density distribution for run m05 for three snapshot times indicated in the title of the individual images. 
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Figure 5. Logarithmic density distribution for run m30 for six snapshot times indicated in the title of the individual images. Our FR-index is above 1.5 for the 
plots at times 4, 14 and 19, and below 1.5 at times 7, 22 and 26. 




Figure 6. Comparison of thermal pressure (red dashed line) against forward ram-pressure (black solid line) on the axis of symmetry for the final snapshot of 
all simulations. 
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M05: div(v) 




50 100 150 

X/L, 



M30: div(v) 




m30: div(v) 




Figure 7. Examples of the divergence of the velocity field, shown is log(| V • v|) for runs M05 (top), M30 (middle) and m30 (bottom). The final snapshot (same 
as in Figure [T) is chosen. Shocks are visible as sharp linear features. Adiabatic expansion (and compression) appears as diffuse structures. The expansion 
region in the initial conical part shows up clearly. 



jet with a significant cocoon, which is confirmed by the simula- 
tion. All the low Mach-number jets have expanded considerably 
beyond L2 by the end of the simulation. Consequently, very weak 
bow shocks, if any at all, are found, except near the tip of the su- 
personically propagating runs m05 and ml 5. In a similar fashion to 
the high Mach-number case, only the simulation with a half open- 
ing angle of 30° develops a stationary terminal shock (Figure |5). 
The predicted position (Lie in Table |2]) is 3.5Li. From the simula- 
tion, we measure a value of l.VLi, and it oscillates between \ALi 
and I.8L1. The reason for the slight disagreement is the pressure 
enhancement downstream of the shock (Figure |6j. 



The cocoon in run m30 has transformed in an interesting way 
(Figure [5]). While for the other simulations, the cocoon is oriented 
backwards, surrounding all of the jet, here it transforms into a 
roughly conical forward flowing structure. Compared with M30, 
in which the cocoon is overpressured by a factor of ^ 100, m30 has 
expanded to more than three times L2 and is in pressure equilib- 
rium with the environment. Hence, if a cocoon had formed around 
the inner jet, it would have become Rayleigh-Taylor unstable. We 
have argued in Section|2]that such a morphological transformation 
should be related to strong entrainment. The re-collimation shock 
is unstable in the high opening angle runs; additionally in run m30 
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M05: integrated compression 



rnOS: integrated compression 
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M15: integrated compression 
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Figure 8. Negative part of the divergence of the velocity field integrated along a path perpendicular to the symmetry axis, i.e. — /^d_y V • v|<o, normalised to 
the respective maximum value. This is a measure of the shock strength. The final snapshot is shown for each simulation. We use ten bins along the X-axis. 
Horizontal error bars are the bin width. Vertical error bars are one standard deviation. The regions of the strongest shocks are where the error bar of the standard 
deviation is largest. 
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the ambient gas is in direct contact with the jet and efficient entrain- 
ment indeed happens at the re-colUmation shock. 

The terminal shock Umit should be determined by the for- 
ward ram pressure becoming equal to the local thermal pressure. 
We confirm this from our simulations in Figure |6] For the runs with 
half opening angle 30°, where the terminal shocks occur before re- 
collimation there is indeed a close match between the two near their 
terminal shock regions (X = X.lLi and X = 11 Li for runs m30 and 
M30, respectively). In run M30, the terminal shock is clearly kept 
at its position by the cocoon pressure. For the other simulations, the 
ram-pressure gets sometimes remarkably close to the thermal pres- 
sure near the position where the re-collimation shock hits the axis. 
It is, however an oblique shock, and no exact match is required. The 
internal jet Mach-number is related to the ratio of the ram pressure 
to the thermal pressure. From Figure |6] the internal Mach-numbers 
are therefore of order ten for runs m05, M05 and Ml 5. For run ml 5, 
it varies between two and six. Even for the high opening angle runs, 
the internal Mach-number downstream of the terminal shock quite 
often exceeds two. 

4.2 Fanaroff- Riley classification 

The simulations we present above show a clearly distinct mor- 
phology. In order to link them to the observed structures of radio 
sources, we have to derive properties of the associated synchrotron 
emission, which is due to non-thermal, relativistic particles and 
magnetic fields. Neither is included in our simulations. For the FR 
classification, we concentrate on low frequency radio emission. The 
latter is however commonly approximated by j p^-^^ where p is 
the thermal pressure (e.g. Saxton et al. 2002). The idea behind this 
is the following: The emitting particles are accelerated in shocks, 
which also raise the pressure. After acceleration, the particles suffer 
adiabatic losses, while being advected away from the shock. Sim- 
ilarly, the pressure in the hydrodynamic simulation declines when 
the plasma re-expands behind a shock wave. The exponent follows 
from the equipartition argument. 

For our low opening angle runs, strong shocks map on to 
pressure enhancements. However, for the high opening angle runs 
(M30, m30), the strongest shock for much of the simulation time 
is the standing terminal shock in the re-collimation region, but the 
pressure is fairly uniform from the post-shock region, downstream 
into the cocoon (Figure [6]). It is therefore more realistic to include 
a measure of the shock strength directly into the emissivity. The 
compression rate in the flow is characterised by V • v, where v is 
the velocity vector. It also traces adiabatic compression and expan- 
sion, but the compression is strongest in shocks. 

Figure |7] shows V • v for runs M05, M30 and m30 for the fi- 
nal snapshots. Besides the initial adiabatic expansion in the conical 
part of the jet, the plots are dominated by linear shock features. 
Clearly, for run M05, the strongest shock is the classical Mach disk 
near the head of the source. Moderately strong shocks are found 
around the Mach disk and somewhat upstream. The shocks are con- 
centrated in the outermost half of the jet. This is in sharp contrast 
to the high opening angle runs M30 and m30. Here, the strongest 
shock is the standing terminal shock in the re-collimation region, 
and many weaker shocks are located in the vicinity. We quantify 
this in Figure [U where we plot the integral of the compression rate 
perpendicular to the axis, — JydyV ■ v| <o (i.e. positive values which 
are not related to shocks are set to zero). We use ten bins along each 
jet in order to avoid overemphasising isolated features, and give the 
bin size and the standard deviation of the compression rate in each 
bin as horizontal and vertical error bars, respectively. All runs with 



half opening angle of 15° or less display the strongest shocks in 
the head region, which is evident from the large error bars there. 
In contrast, for the high opening angle runs the standing terminal 
shock in the re-collimation area is the strongest feature, also, when 
averages over larger areas are considered, as we have done here. 
Correspondingly, the compression rate reaches the highest values 
near that site, in the innermost part of the jet. Thus, the shock struc- 
ture in high opening angle runs is clearly distinct from the one in 
low opening angle simulations. Yet, the maximum of the pressure 
is always close to the tip of the jet, also for runs M30 and m30. 
We exemplify this for the final snapshot of run m30 in Figure |9] 
where we compare emission maps using p^-^ and V • v| <o as respec- 
tive emission tracers. Clearly, the assumption that pressure maxima 
trace strong shocks is not generally correct for our high opening 
angle runs, and we need another proxy for the emissivity. 

Shocks are most directly traced by the compression rate V • v. 
Realistically, the particles would suffer losses while streaming 
away from the shocks, which smooths out the emission over a cer- 
tain length behind the shock. We cannot hope to model the de- 
tailed appearance of these emission features with our simulations in 
any useful way. However, the total emission initiated by any given 
shock should be related to the strength of the shock and its area. 

For connecting our results to observations we therefore be- 
lieve it is appropriate to take the integrated compression rate de- 
fined above (— J^dj V ■ v|<o) as a measure of the emissivity. At the 
same time we have reduced the spatial resolution so that the emis- 
sion due to each shock is assigned to a larger volume. 

Also in observations the radio emission appears knotty. A 
given radio knot is most likely retaled to one particular shock or 
shock complex. Additionally, the FR classificat ion has originally 
been established using low resolution radio maps (iFanaroff & RilevI 
Il974l) . It is therefore appropriate that we also smooth the simulation 
results to a comparable resolution. 

With this choice for the emissivity, we produced pseudo-radio 
maps for all the 2949 snapshots of our six runs, and smoothed it to 
20 resolution elements along the axis. We defined the source size by 
the extent of the pseudo-radio emission along the X-axis, up to the 
point where the emission drops below five per cent of the maximum 
value (xsize)- We also calculate the position of the brightest feature 
^bright- Finally, we define the FR-index as: 

FR = 2xbright Asize + 1/2. (10) 

With this definition, 0.5 <FR< 1.5 for an FR I radio source, and 
1.5 <FR< 2.5 for an FR II radio source. 

We show FR against simulation time in Figure [TO] For 
runs M05 and Ml 5, the average FR is almost always clearly above 
1.7, indicating an FR II source, as expected. The values below 1.5 
are mostly due to shocks near the inner boundary, which sometimes 
get quite strong, especially because of the reflecting boundary con- 
ditions we use. This artifact of the procedure does however not de- 
stroy the general trend. We also find a clear FR II morphology for 
the heavy low opening angle jets m05 and ml 5. Run m30 is an 
FR II source until it reaches its terminal shock limit Li\y around 
t = 4. From then onwards, it changes its character to an FR I radio 
source. For run M30, the FR index is on average above 1.5 for times 
earlier then t = 0.1, and below afterwards - but only marginally so. 
So, it makes principally the same transition as run m30, but the 
transition is much less marked. In this case the jet is quite un- 
derdense, forming a prominent cocoon, which protects it against 
entrainment of the ambient medium; it therefore transports more 
energy to its head region than the m30 jet. Consequently, the head 
region also contains strong shocks, and the emission maximum is 
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often located there. We expect that the FR II branch of the diagram, 
which is already weaker than the FR I branch after t = 0.7, would 
die out completely, as in run m30, if there were stronger entrain- 
ment. This could happen either when the cocoon approaches the 
pressure of the environment, or if it had a sufficient source of gas, 
for example a gas disk in its equatorial plane. 



5 DISCUSSION 

We have demonstrated both analytically and by hydrodynamics 
simulation that for jets collimated by the external pressure, the 
large-scale morphology developed as a result of the interaction of 
the jet with its surroundings near its collimation region. The ratios 
of three critical length scales to one another control the behaviour. 
Only two of them are independent, and they are set by the jet's 
opening angle and external Mach-number. 

We find that light jets with cocoons form for high external 
Mach-numbers and moderate to large opening angles. It is well 
known that the majority of extragalactic radio sources have ex- 
tended radio lobes and the jets are therefore very light compared 
to the external medium. 

We predict that all lobed sources go through a phase in which 
hot spots occur near the tip of the lobe. If the jets collimate before 
their ram-pressure falls below the thermal pressure of the environ- 
ment, then the jet retains FR II morphology. In our combined an- 
alytical / simulation analysis, this depends on the initial opening 
angle of the jet with a transition for 6 = 24°. Our simulations with 
5° and 15° confirm this scenario, and consistently develop FR II 
morphology. 

The jets with 30° run out of thrust, before the re-collimation 
shock reaches the axis. This result does not depend on the external 
Mach-number, which differs by a factor of 100 between our two 
30° runs. 

Run m30 has high opening angle and low external Mach num- 
ber. This results in a collimated jet density that is only slightly be- 
low the ambient density. Consequently, the cocoon is not as pro- 
nounced as in the higher external Mach number runs, and is not 
able to shield the jet from the ambient medium. Instabilities at the 
re-collimation shock, which are pronounced at high opening angle, 
entrain ambient gas directly in the jet in this run. This slows the jet 
down, which produces a conical turbulent flow downstream of the 
re-collimation shock. Because of the large opening angle, the ram 
pressure of the jet is insufficient to propagate the terminal shock 
through the re-collimation shock. In effect, the jet shows FR II mor- 
phology until the terminal shock reaches its limiting location, and 
appears as an FR I source afterwards. 

Run M30 has a very similar shock structure to that of run m30. 
We also classify it on average as FR I at late times. This suggests 
that the transition from class II to class I is independent of the jet ve- 
locity, which differs by two orders of magnitude between the runs, 
but only depends on the initial opening angle. 

The cocoon of run M30 differs from that of run m30: it has 
lower density and at the end of the simulation it still effectively 
prevents entrainment of ambient gas. As argued above, this also 
leads frequently (but still in the minority of snapshots) to dominant 
shocks in the head region. M30's cocoon could transform to a mor- 
phology similar to that of run m30 if it entrained sufficient amounts 
of dense ambient gas. The cocoon itself is stabilised against entrain- 
ment of a mbient gas due to Rayleigh-Taylor instabilities by its de- 
celeration (I Alexande J2002h . The deceleration up to a source size of 
L2 can be approximated by a power law (e.g. iKaiser & AlexandeJ 



Il997h . This power law breaks down when the source ap proaches 
L2, and the sideways expansion stalls (iGaibler et al.l200 3 show this 
for FR II jet simulations). The cocoon is then Rayleigh-Taylor un- 
stable and will be refilled with ambient gas, which will also lead to 
entrainment into the main beam (i.e. the forward flow downstream 
of the re-collimation shock). The cocoon will transform accord- 
ingly and this should result in a source structure similar to the m30 
run. Our M30 run is however a factor of 100 away from this tran- 
sition point, and for the present stud y we did not have the compu- 
tational resources to simulate this far. iKaiser & BestI (l2007h invoke 
this expansion to L2 and the subsequent entrainment as the sole 
mechanism to account for the FR dichotomy. 

Observations of the FR I source Centaurus A (e.g. 
iHardcastle et al.ll2007b suggest another source of entrainment. If 
the re-collimation shock is located inside a dense gas disk, that gas 
disk would be difficult to push aside, and the dense gas might come 
close to t he beam jNeu maver et al. 2007), especially if the latter 
precesses (iMorganti et al.ll 19991). Entr ained gas can be resolved in 
this source fe.g. lTingav & Lendl2009h . This confirms that entrain- 
ment might be an important factor for a radio source to show FR I 
morphology. Centaurus A shows^vidence for a high Mach num- 
ber bow shock (.Croston et al.ii2009i) , and is therefore still far from 
pressure e quilibrium with its s urroundings. A mechanism as pro- 
posed by (IKaiser & Bestll2007h . which relies entirely on the radio 
source to be close to pressure equilibrium with its surroundings, 
therefore seems to have difficulties to explain the observations of 
Centaurus A. Stellar winds within the jet dKom issarov 1994) are 
another possible mechanism by which dense gas is brought to the 
jet plasma. 

While entrainment has been recognised in many studies as be- 
ing important for the emergence of an FR I morphology, we find 
there is another, necessary condition - a large opening angle for the 
jet. Regarding observed opening angles, proba bly the best s tudied 
cases are M87 (FR I) and Cygnus A (FR II). ljunor et all (Il999l) 
present an observational study of the inner jet in M87. The jet 
starts out at an apparent half opening angle of about 30° at about 
0.04 pc distance from the core. Its opening angle then reduces 
gradually to about 7° about 14 pc distance from the core. In our 
simulations, collimation also happens gradually, over a comp arable 
length scale ratio of about two orders of magnitude. Gracia ^taD 
(l2009h present a detailed analysis of the inner jet of M87, using 
a magnetohydrodynamics-based fit to the observational data. The 
available data constrains their models well. Among 2600 MHD 
models, they find several that fit the data well, apart from the in- 
nermost point with the largest opening angle. Discarding the inner- 
most data point yields a second innermost opening angle of only 
about 17°. While this might first appear to contradict our results 
(see below), the other features of their best-fitting models are in 
good agreement. They find that the limb brightening of the jet is 
due to enhanced intrinsic emission, rather than a beaming effect 
due to velocity structure. Further, they identify the bright feature 
'knot A with a re-collimation shock. If the jet density is related 
to the emissivity, limb brightening before the re-collimation shock 
reaches the axis is predicted by our simulations. It is a natural con- 
sequence of the transition from a co nical to a pressure collimat ed 
jet, which has already been noted bv lKomissarov & Fallel(ll998h . 

The opening angles derived from observations cannot be di- 
rectly compared to the initial opening angles we impose in our sim- 
ulations or use in our analytic work. This is because the jet base is 
both hard to resolve and might not display bright features over the 
full jet width. The first reflection point of the re-collimation shock 
on the axis, or for FR I sources the flaring point will be locations of 



12 M. Krause et al. 



nn30: ray-traced emission: p' 



mSO: ray-traced emission: div(v) 




20 1 



10 



-10 



-20 I 



10 20 30 40 



IP 20 30 



40 



Figure 9. Line of sight (perpendicular to the jet axis) integrated, 3D transformed emissivity for different emission tracers. 20 resolution elements are used for 
the X-axis, p^-^ (left) badly reproduces strong shocks. | V • v| (right) highlights shocks much better. We use the latter for comparison to observations. See text 
for details. 



bright emission, where an opening angle could easily be measured. 
Using this point, we obtain measured half opening angles of 5° 
for our simulations of 30° initi al opening angle. The half opening 
angle at knot A in M87 is 3° dJunor et all 1 19991) . The agreement 
with our results is reasonably good for M87, given that the ob- 
servational value is a lower limit, since the full width at quarter 
maximum is used, and in addition, our simulations do not take into 
account magnetic and relativistic effects. 

Measuring the opening angle at the first reflection point of 
the re-collimation sho ck yields a half ope n ing an gle of ~ 1° for 
runs M05 and M15. iKomissarov & Fallel (1 19981) present similar 
simulations to our low opening angle cases and compare then in 
detail to the FR II radio source Cygnus A. From Cygnus As mea- 
sured cocoon aspect ratio, they predict an initial half opening angle 
between 11.5° and 15.25°. They identify knot 3 at 8 arcsec (about 
6.4 kpc) distance from the core as the first reflection point of the 
re-collimation shock. As in the case of M87, they mention the limb 
brightening to support this identification. They quote an opening 
angle of 0.8°, measured at 20 arcsec from the core, pointing out that 
the measurement is inconsistent with their predicted initial opening 
angle. However, when doing the measurement in a similar way for 
the simulation and for the observation, the opening angle on the 
scale of the re-collimation shock is actually very similar. 

The most important parameter that correlates with the ob- 
served jet morphology is the jet power. We find that only jets with a 
large initial opening angle have the potential to form an FR I source. 
This s uggests an anti-correlati on of the opening angle with the jet 
power. jPushkarev et al.l (l2009l) present a sample of parsec scale ra- 
dio sources, all having opening angles of order a few degrees. In- 
terestingly, they find an anti-correlation between the opening angle 
and the L orentz-factor, which is co nsistent with the jet acceleration 
model of iKomissarov et all (l2007l) . Therefore, there might indeed 
be a correlation of a jet's opening angle with its power. We are not 



aware of a systematic study of the correlation of the jet opening 
angle with FR class. 

We identify the quasi- stationary terminal shock in the high 
opening angle runs with the flarin g point in FR I jets. This is in line 
with other theoretical w ork (e.g. lBicknelliil984 ; jPemcho & Marti 
i2007l) . Observationally, iHardcastle et al.i (l2005i) and iJetha et al. 
([20051) find that the flaring point is located in a region of flatten- 
ing ambient pressure profile - a result that can be understood in 
the sense that in such a region it is much more likely that the jet's 
ram-pressure will become equal to the ambient pressure. 

If the identification of the flaring point with Lie is correct, 
the location of the flaring point would constrain Li . Together with 
information on the jet speed, the external density and pressure, this 
would provide the possibility of measuring the power of FR I jets. 

We do not expect that the general features we discuss here 
would change if we had do ne the simulations in three dimensions. 
iKomissarov & Fall^ (il998i) have done simulations of low open- 
ing angle axisymmetric FR II jets for both the relativistic and the 
non-relativistic case, with very comparable results. The critical an- 
gle for the morphological transition depends on the angle that the 
re-collimation shock makes with the axis. In this respect, our re- 
sults for 5° and 15° opening angle are quite similar to their 5° 
and 20° cases. Also their relativistic simulations show quite simi- 
lar behavio ur. The critical angle of 24° we derive is consistent with 
the fact that iKomissarov & Fall9 (il998i) still find FR II structure at 
20° for their relativistic and non-relativistic cases alike. While the 
value of the critical angle might depend on where the re-collimation 
shock exactly meets the axis, which may well depend on details and 
also possibly on the magnetic field, we expect the general result to 
be robust. 




Figure 10. Fanaroff-Riley class over time for all simulations. The solid line with error bars denote the mean value and standard deviation for the given interval. 
Additionally, all snapshots are represented by filled circles. The red line denotes the border between FR I and FR II. See text for more details. 
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6 CONCLUSIONS 

We have shown that the large-scale morphology of pressure colli- 
mated jets is determined by the processes near the region where the 
jet collimates due to the ambient pressure. We find that the main 
requirement to form a very underdense jet, and therefore a wide 
radio lobe is a high external Mach number. We find a clearly dis- 
tinct shock structure for jets with initial half opening angles below 
and above about 24°, respectively. Jets with lower opening angles 
always have FR II morphology, and we classify them accordingly 
by an appropriate FR-index based on the shock structure. Higher 
opening angle jets are found to transition from class II to class I. 
In the latter case, entrainment is identified as an important factor in 
the jet morphology. The initial opening angle near the jet base is 
much larger than the observed one. Predicting the observed open- 
ing angle from our simulations results in values which agree well 
with observations of M87 and Cygnus A. 
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